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THE DISTRIBUTION OF CURRENT ON CYLINDRICAL RADIATING DIPOLES 



SUMMARY 

The distribution of current on a dipole comprising cylindrical conductors 
is approximately sinusoidal, but in certain problems a more accurate solution is 
required, A number of theoretical attempts have been made to obtain a better 
approximation and are in reasonable agreement. These methods are, however, laborious 
to apply in practice and are also limited in application. 

In this report a simple method is described which enables a current dis- 
tribution more accurate than the sinusoidal approximation to be determined. The 
method is based on the perturbation principle and besides being applicable to a wide 
variety of cases, is in good agreement with more rigorous theoretical approaches and 
*dth measurements. 



1, IN1R0DUCTI0N 

A radiating dipole comprising cylindrical conductors behaves approximately 
as an open-circuited uniform lossless balanced transmission line as regards the 
distribution of current along its length. The current is thus distributed approxi- 
mately sinusoidally and the radiation pattern can be calculated on this basis. For 
many applications, the assumption of a sinusoidally distributed current results in a 
calculated radiation pattern which is sufficiently accurate for practical purposes. 
In certain cases, however, an assumption of a sinusoidally distributed current can 
lead to appreciable error in the radiation pattern and a better approximation is 
needed. The calculation of a more accurate distribution is, however, a matter of 
some difficulty since it is not possible to state the field conditions at the tips of 
the dipole precisely. 

It is convenient to assume that the true dipole current is the sum of a 
sinusoidal primary current (a standing wave pattern with a velocity of propagation 
equal to that of light) and a correcting tern which may be called the secondary 
current. The secondary current may be resolved into two components, one in phase 

and one in quadrature with the primary current. The primary current and the in-phase 
component of the secondary current may be regarded as a modified primary current and 
result in a standing wave pattern viiich is substantially sinusoidal but with avelocity 
of propagation less than that of light. The quadrature component of the secondary 
current can conveniently be called the "feed" current^ since it is in phase with the 
driving voltage and is associated with the radiation of power. 

Theoretical attempts have been made " to obtain an approximation to the 
secondary current and are in qualitative agreement on the manner in which the true 
current distribution varies from a sinusoidal distribution. However, these methods 
are laborious to apply to practical cases and difficulties arise in the case of a 
dipole vAiich is mutually coupled to other aerials. 



In this report a method based on the perturbation principle is described 
viiich enables a current distribution more accurate than the sinusoidal distribution 
to be determined. The method is simple to apply and gives a result for the "feed" 
current viiich is in good agreement with that derived from more rigorous approaches. 
In addition, the solution of cases in which the aerial is mutually coupled to other 
aerials is straightforward. The extension of the method to predict the in-phase 
component of the secondary current is not simple, but this is not thought to be 
a disadvantage. The theoretical analyses^"® are not in amy case in good agreement 
as regards prediction of the in-phase component and reliable measured figures of the 
reduction of velocity of propagation are available. ^ 



2. BASIS OP METHOD 

2. 1. The Perturbation Principle 

It has been pointed out* that the compensation theorem for elefctrical 
networks, though capable of exact statement, is principally useful as a perturbation 
method. Such a method enables the current flowing in any conductor in a linear 
reciprocal network to be determined from the change of input impedance vAiich occurs 
vAien a small perturbing impedance is placed in the path of the current. In many 
cases it is much simpler to determine the change of input impedance than to calculate 
the current directly and the method is thus of considerable interest. 

If I is the current flowing in a conductor in a network and an impedance 
82 is inserted in series with the conductor, the input impedance will change by an 
amount SB (say). It can be shown that: 

I 

Bz = —(I + 81) Bz (1) 

-To 

where Sj is the change of current resulting from the insertion of the perturbing 
impedance 82 and Iq is the current applied to the input terminals of the network. 
Equation (1) is exact but provided that 82 is made sufficiently small it can be 
simplified to:* 



1= Bz 

lo S2 

Thus the current flowing in any part of a network can be determined from 
the change of input impedance due to a perturbing impedance. For the purposes of 
analysis it is usually convenient to assume that Iq is unity. 

2.2. Perturbation of a Cylindrical Dipole 

The perturbation method described in Section 2.1 gives a means of obtaining 
an approximation to the secondary current on a cylindrical dipole; for convenience, 
however, the electrically identical case of a unipole above a perfectly conducting 
ground plane will be considered. If the unipole is broken at any point and a small 
impedance is inserted in series (for convenience a small inductance will be assumed), 
then the current at the perturbed point can be determined from the change of input 



impedance. As an initial approximation the current distribution is assumed to be 
sinusoidal and this enables an approximation to the input impedance to be readily 
calculated. This input impedance will be: 

Z = R -t jX (3) 

where R is the radiation resistance (determined by integration of power over a large 
sphere) referred to the input terminals and /is the input reactance -Zq cot Zrrl/K 
the theoretical value for a lossless transmission line. Zq is the characteristic 
impedance of the unipole* and I the length. 

When the perturbing impedance is inserted the current distribution is 
changed. If this change of current distribution is calculated on the assumption 
that the unipole behaves as a lossless transmission line the change of input impedance 
will be: 

S^ = hR + jhX (4) 

where Si? is the change of radiation resistance due to the change of current distri- 
bution and 8/ is the change of input reactance of a loss-free transmission line. 
■fiius from the change of impedance 8^ given by equation (4), an approximation to the 
current in the unipole at the point of perturbation can be derived from equation (2). 
This current will differ from the original sinusoidal assumption since account has 
now been taken of the radiated field by including the change of radiation resistance. 

The error in the determination of the secondary current will clearly depend 
on the accuracy of the initial assumption of sinusoidal current distribution. Thus 
considerable error will result if the perturbation method is applied to regions of 
the unipole in which the sinusoidal primary current is small; for example, at or 
near the drive point of a unipole K/2 high. At the top of a unipole both the primary 
and secondary currents must of necessity tend to zero and the error in the secondary 
current determined is difficult to assess. The current at the top of a unipole 
camnot, however, be determined with accuracy by any method due to the difficulties 
involved in assessing the field in the neighbourhood of the tip. 



3. CDRREaTT DISTRIBUTION OP UNIPOLES OF HBIfflT UP TO k/A 

The case considered is shown in Pig. 1(a), the unipole having a characteris- 
tic impedemce of Zq. The height is h and Znh/K - cf^Q vAiere k is the free space wave- 
length. Pig. Kb) shows the equivalent lossless transmission line from which the 
preliminary assumptions of current distribution before and after perturbation are 
determined; it is assumed that unit current is injected at the base from a constant 
current generator. 

'nj&llke a Qon-radiat lag transmission line a unipole cannot hare a characteristic iKpedance since 
it is not uniform. HoveTer^ for practical purposes, the unipole bebares approximately as a 
uniform transmission line and a parameter corresponding to the characteristic impedance can be 
used. Of sOTeral formulae, the one In best accord ¥lth experimental results is due to Q.T.O. 
Howe^'' and is Zq — 60(logj ( i/'*' — 1) ohms where I is the length andcctbe radius of the cylinder. 
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The unipole is perturbed at a distance of 6 radians from the top by a small 
series impedance jtlZq and on the assumption of a lossless transmission line the change 
of input reactance (equation (4)) can be shown to be: 



jhX = jZd^ 



sin ^00 



(5) 



Before the perturbation is inserted the current at a point <p radians from the top of 
the unipole will be: 



lp(0) = 



sm ' 



sm (pQ 



(63 



where the subscript P denotes the sinusoidal primary current. After insertion of 
the perturbation the current at a point (p radians from the top of the unipole will be: 



_ sin(g + A) 
^ sin((J^ + A)sin^ 



sin(0 + A) 
^ sln{4>o + A) 



3in<^ 



0< 0< 19 



(7) 



vi.ere A= n sin d the effective change of length of the unipole in radians. 



If the height of the unipole is sufficiently small, the vertical radiation 
pattern (v.r.p.) will be unchanged when the perturbation is added. The horizontal 
radiation pattern (h.r.p.) will be circular irrespective of the perturbation and the 
change in radiation resistance can thus be determined from the total current moments 
of the unperturbed and perturbed aerials. These current moments will be Mi for the 
unperturbed unipole and M^^ for the perturbed unipole and will be determined by inte- 
gration of equations (6) and (7) giving: 

*<^o 
^ ' ^c 



M^ = 




sin<i (ii = tan -i-2 



M n — 



sin( 9 + A) 
siii(c/)o + A)sin6' 



■e 



sincf) dcfi + 



sx-nicpQ + A) 



f4>o 

sin(c^ + A) 



•/O r- 



(8) 



tan^ 

2 



1 + 



A(tan^ - tan4) 

1 - cos 1^0 



Thus if the radiation resistance of the unperturbed unipole (measured as 
a series component at the drive point) is R the radiation resistance of a perturbed 
unipole is given by: 



(9) 





M% 




ff' 





^u 


R 


Ml 


H 



Substituting for Mi, M2 and A from equations (7) and C8) the change of 
input resistance due to the perturbation is given by hit = R' - R and is: 



n sin^i9(tan^ - tan-|) 

" ' 71;^ -'' 



(10) 



Equations (5) and (10) give the change of input impedance due to the pertur- 
bation, and from equation (2) the square of the current at the point of perturbation 
is given by dividing the change of impedance by the value of the perturbation (juZq); 
thus: 

2 S^ 8/? + jSx 2 j 2 , ^D 2 <5^o J, a \ 

I (6) = —— = = sin 9 <eosec <Po - j —^ cosec (tan?p - tan -5-) > 

hz jnZo (_ Zq 2 ^ '^ j 

where the ahsence of the subscript P viiich was used in equation (7) denotes that the 
current is now the total current comprising both the primary and secondary terms. 
It should be noted that the primary currents of equations (6) and (7) and the total 
current of equation (11) appear as functions of two different variables. The 
primary current which is a function of 4> must be integrated over the length of the 
unipole to find the total current at the perturbed point 9 radians from the top of the 
aerial. 

In most practical cases, the imaginary term in equation (11) has a modulus 
small compared with unity and the current may be determined by approximating to the 
square root giving: 



1(9) = s±n9 Jcosec</)o 



j— (1 - tan-^cot 



-tnnt.^O) 

2 



(12) 



where the real term is the sinusoidal "primary" current and the imaginary term is the 
quadrature "feed" current. The modified current distribution must account for the 
power radiated from the unipole; the sinusoidal primary current fails to do this 
since it results in a purely reactive input impedance. The perturbation method 
assumes that the unipole behaves as a uniform transmission line and the voltage at 
a point 9 from the top is therefore 7(9) - -jZq 'dl(9)/W. Applying this relationship 



to equation (18), the voltage at = 0o (i-e. at the drive point) is R ^- j ZQCotcpQ. 
Since an input current of unity has been assumed the input impedance is R^-jZocoi/j^ 
and the correct radiated power is therefore drawn from the generator. 

The relative and not the absolute amplitudes of the two terms corresponding 
to the dipole current are of interest and it is more convenient to write the sinu- 
soidal terra without a multiplying constant. Thus equation (18) becomes: 

/ R tanf \ 

i^ [6) = sini9 [1 - j -L a - — T-) sing^b I (la) 

I 2o tan'^ j 

vAiere r denotes the relative current. 

As has already been pointed out, equation (13) is based on the assumption 
that vdien the perturbation is inserted in the dipole, the change in the shape of the 
v.r.p. can be neglected. In the Appendix it is shown that this assumption is valid 
for unipole heights of up to \/A and in this case cpQ = tt/S and equation (13) becomes: 

lid) ■= sin^ 1 - j_£ (1 - tan-|) (14) 

Another case of interest is that of a very short unipole with <^o « n/Z in 
viiich case equation (13) becomes: 



I, m 



= ^ |i - jj!* (<^ - e\ (15) 



In applying equations (13), (14) and (15) in practical cases, values of ^„, 
the radiation resistance of the unipole referred to the drive point are necessary. 
Figures for the radiation resistance based on a sinusoidal current distribution for 
dipoles of arbitrary length are readily available. ^^ 



4. CURRENT DISTRIBUTION IN THE PRESENCE OF EXTERNAL FIELDS 

In Section 3 the current distribution on an isolated unipole was analysed. 
The case now considered is that of a unipole in the presence of a second aerial of 
arbitrary type illustrated diagraimiatically in Fig. 8. Both aerials are driven from 
a common source delivering unit current, a current of I-i being impressed at the 
unipole input terminals and I2 at the terminals of the second aerial. An equivalent 
circuit of the arrangement can be drawn as in Fig. 3 where Z^ and ^2 a^re the self 
impedances of the unipole and second aerial respectively and Z^ is the mutual impe- 
dance.* 

It is required to find the current on the unipole and a series perturbation 
of jn^o is therefore inserted as in the analysis of Section 3. The insertion of 
the perturbation will cause both the self and mutual impedances of the unipole to 

Z^ and Z r^ are the Impedances of each aerial with tbe terminals of the other aerial open circuited 
and Z is tbe open circuit voltage at the terminals of one aerial with 1 ampere Injected into 
the other. Tbe second aerial is assumed to be such that open circuiting Its terminals will 
result In Its effect on the unipole being negligible. 
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change; suppose these changes are S^^ and hZ^ respectively. The equivalent circuit 
then becomes as in Pig. 4 and from the theorem of equation (8) the change of input 
impedance seen by the common source will be: 



IZ = llS^i + Slilg S^, 



(16) 



Also from the theorem of equation (2) the square of the current on the 
unipole at the point of perturbation must be: 






(17) 



where k is the ratio of the currents in the aerials, i.e. Ig/Ii. 

^y definition S.2i, the change of self impedance of the unipole, is indepen- 
dent of the second aerial and is given by equations (5) and (10), The change in 
Z^ will be given approximately by: 



IZ^ = *n^. 

Equation (18) assume^ that 
the perturbation changes the amplitude 
of the field radiated by the unipole 
but not its v.r.p. This is the same 
assumption which has already been made 
in computing the change of input 
resistance of an isolated unipole and 
is discussed in the Appendix. 

Prom equations (5), (10), 
(17) and (18) the square of the current 
distribution on the unipole is thus 
identical with that on an isolated 
unipole in which R- is replaced by 
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(tan^ - tan—) 



(18) 
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(S^ + kZ. ) 



Provided that (R„ + 'kZ ) is sufficiently small compared with 2n, the 
square root can be approximated and the expression for the current is that given in 
equation (13) with (R. + k^_ ) substituted for if„. The mutual impedance will in 

ft n 

general be complex and thus both the in— phase and quadrature components of the total 
aerial current will be modified. In the case of the \/4 unipole the current distri- 
bution for a unit input current to the unipole will be: 



1(6) = sin6 



(R^ + k^.)(l - tanS.) 



(19) 



5. CURRENT DISTRIBUTION ON A HALF-WAVE UNIPOLE 

The case of a \/2 unipole is of considerable interest since it forms the 
basis of the m.f. "antifading" type of transmitting aerial. The perturbatibn method 
outlined in this report cannot be applied directly to this case since the sinusoidal 
first approximation of current is zero at the drive point and a large error would 
occur in the calculation of the "feed" current. The error occurs because the 
necessary condition that the secondary current must be small compared with the sinu- 
soidal primary current, is not satisfied. The current distribution already deter- 
mined for a A./4 unipole can, however, be used to obtain the feed current distribution 
for a K/2 unipole. 

Fig. 5 shows the case considered and the initial sinusoidal approximation 
of current; the "loop" current is assumed to be unity. The power radiated can be 
determined by integration of the distant field and since the loop current is J amp 
it will be numerically equal to R , the radiation resistance referred to the current 
loop. On the initial assumption that the aerial behaves as a lossless transmission 
line, the input voltage must be ^q volts and since the input power is equal to ^f, the 
input current must be R /Zq. It should be noted that the input current must be 
advanced in plane by tt/2 radians relative to the loop current and an approximation 
to the true current, taking into account radiation, is therefore known at two points 
on the aerial, the mid point (the current loop) Eind the drive point; this is illus- 
trated in Fig. 6. The current distribution will not be changed if an infinite 
impedance generator giving unit current is inserted at the current loop as illustrated 
in Fig. 7. The current on the unipole can thus be found by the superposition of 
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the currents obtaining in the two 

arrangements of Pigs. 8(a) and (b). 

In computing the current distribution 

of arrangement 8(a) equation (19) is 

used, the values of Z, being determined 

from the formulae for mutual impedance 

between half-wave dipoles given by 

Carter* (the value of mutual impedance 

thus derived is 2^ = 26-4 + j20-Z 

ohms). In computing the current 

distribution for Fig. 8(b) the current 

on the section marked A is ignored 

since it will be small. The resulting 

current distribution for a A./S unipole 

of characteristic impedance of 250 ohms is shown in Pig. 9 together with the result of 

a calculation made by G.D. Monteath based on BOhm's theory; the agreement is seen to 

be good. Pigs. 10 and 11 show comparisons of measured and calculated results for 

the amplitude and phase of the current on aA./2 unipole and once again the agreement is 

seen to be good. 
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6. RESONANT LENGTH OP DIPOLES AND UNIPOLES 

It has been mentioned in the Introduction that cylindrical unipole and 
dipole aerials resonate (i.e. have zero input reactance or susceptance) for lengths 
less than multiples of A./4 in free space and that this effect is due to the in-phase 
component of the secondary current. The perturbation method of determining the feed 
current depends on a sinusoidal first approximation to the current distribution; this 
need not necessarily assume a velocity of propagation equal to that of light. The 
initial sinusoidal current distribution may be conceived to correspond to a reduced 
velocity of propagation and this will clearly result in a more accurate prediction 
of the total current; the calculated curves in Pigs. 10 and 11 have been derived in 
this way. The perturbation method does not seem to lead readily to an estimate of 



-r 



MUTUAL IMPEDANCE, 
BETWEEN OIPOLE AND 
IMAGE EQUALS Z. 






\ 

T 



, CONSTANT CURRENT 
I AMP GENERATOR 



X 



/ /// /y 



( I', IMAGE 



///////// 






CONSTANT CURRENT 
GENERATOR_ 



NJECTING /^AMPS. 
/////// ' 



it) 



(a) 



Fig. 8 - Superposition of generators for hal f- wave length unipole 



10 



10 



0-8 



^ 0-6 



0-4 



d 02 



-0-2 









/'' 




^\ 


Ni 








/ 


/ 

/ 

• 








\ 

\ 






/ 












\ 




/ 


f 












\ 


/ 




^ 


r^~v^^ 










\ 
\ 










**^^ 


^^^^^=^ 


.=^ 




^ 



I80" I60" 140" I20° lOO" 80" 60" 40" 20" O" 



t 



Calculation by 

Perturbation method 

Calculation based 

on Bohm' s theory 

Primary current 



GENERATOR 



* <- 



TOP OF UNIPOLE 



Fig. 9 - Comparison of calculated current on half-wavelength unipole 

reduction of velocity although a figure could be deduced by consideration of a new 
approximate transmission line including the effect of "feed" current. However 
predictions of the reduction of velocity along the aerial given by other theories 
do not agree with measured results and measured figures^ have been used in the cal- 
culations for Pigs. 10 and 11. 
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CONCLUSIONS 



A simple method is given of obtaining a correction to the first order 
assumption of a sinusoidally distributed current on a cylindrical dipole. The 
correction term derived is in quadrature with the sinusoidal current and corresponds 
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to the "feed" current, i.e. the current component associated with radiation of power. 
The method is applicable to unipoles of any length up to K/A (or to the corresponding 
balanced dipoles) and can be extended by superposition to longer aerials. The 
method also enables the current distribution to be assessed in the presence of other 
aerials. 
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APPENDIX 
The Power Gain of a Perturbed Aerial 

In this report, the change in radiation resistance due to perturbation of 
a unipole by a small series reactance is found by assuming that the power gain of the 
unipole is unchanged, i.e. that the radiation pattern is unchanged by the presence of 
the perturbation. This assumption is an approximation and the power gain will in 
fact be modified slightly but the error will be negligible for very short unipoles 
where the v.r.p. is virtually independent of the current distribution. The error 
increases as the unipole is made longer and the worst case considered is that of a 
A./4 unipole. In this Appendix it is shown that even with a \/4 unipole the error 
in the change of radiation resistance is less than 3$ and results in a negligible 
error of the calculated value of the current; the error for lower unipole heights 
will be correspondingly smaller. 

On the initial assumption of a lossless transmission line the current 
distribution on a perturbed ^/4 unipole is given by equation (7) which can be re- 
written as: 

I{4)) = A sln4> + B sin{4> - 9) fSO) 

vAiere -4 = 1 + n sin^ cos^ 

and B = for ^ 4> 4 6 

5 = -n sins' for 6 4, <P 4 ^/2 

Thus the current consists of two sinusoidal terms, one extending along the whole 
length of the unipole and the other extending from the point of perturbation back to 
the input terminals. Because the value of the perturbation jnZo is small compared 
with Zq we have n«l and the second term is much smaller than the first. 

The current on an unperturbed unipole is sin0 and with a suitable choice 
of units the radiated power is: 

I cos ^(77/2 sini/' ) 

Pi = ; d\p = 0-6093 (21) 

I cos \p 

Jo 

where ^ is the angle measured from the horizontal and the solution of the integral 
is obtained from tabulated values of the cosine integral. 

With the same choice of units, the radiated power when the unipole is 
perturbed is: 

^2 = (1 + 2n sinfi" cos6') 



cos {77/2 sinip) 

■ dxfj 



cos 



^ 



- 2a sin^Cl - sin5) 

^o 



^o 
2 

cos(T7-/2 sin0) . F(^) di-p 



^/2 f22' 
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viiere P(\p) is the v.r.p. of the current distribution term B sin(!^ - 6) inequation (20) 
normalized to unity in the horizontal plane. Fi'-p) must lie between the v.r.p.s of 
a A./4 unipole and a very short unipole, i.e. between cos (7t/8 sini//)/cos0 and cos0. 
Hence the value of the second integral in equation (22) must lie between 

cos (77/2 sin (// ) 



cos i/i 



dip 



aad 

COS (77/2 sini// )cos i//.d 0, 

'o 

that is between 0-6093 and 0*6366. The solution of the first of these two integrals 
was given in equation (21) and the second is readily integrable directly. 

If the second integral in equation (22) is assumed to be equal to 0*6093 
the power Pg is that derived on the assumption of unchanged gain when the unipole is 
perturbed. The error in assessing the radiated power is therefore given by: 

Qdj] sin9 (1 - sinj^) (23) 

viiere 17 is a number vftiich is equal to or less than 0*6366 - 0*6093 = 0*0273. 

Prom equation (10) with (po — tt/2 the computed change in radiated power due 
to the presence of the perturbation is (on the assumption of unchanged gain): 

2n sin 6(1 - tan 6/2)Pi (24) 

and hence the percentage error in the change in radiated power is: 

1007] (1 - sin6) 



0*6093 ' [/2 sin ((9 + 77/4) - l] 



(25) 



Equation (25) is zero when 9 = tt/2 and also when 6=0 (because T) = yflien 6=0] 
and the maximum error may thus be expected in the region of 6 = tt/A- If the upper 
limit of 0*0273 is substituted for r], and 6 = 77/4, equation (25) shows that the 
maximum error in the assessment of radiated power is less than; 

100 . 0*0273 (1 - 0*7071) 

. — ■ = 3*25f 

0*6093 (1*4142 - 1) 

From equations (11) and (12) the error in the secondary current found by the pertur- 
bation method is therefore less than Isi viiich is negligible. 
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